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all  dimensions  are  in  mm 


Figure  2:  Crosssection  of  test  models. 
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alt 

(b)  t/R.  =  1.0 

•  Results  from  BEM  showing  “dish”  shaped  region  for  normalized  SIF  for  a 
semi-elliptic  surface  crack  in  a  Dressurized  thick  walled  cylinder. 


•  Ff  -  Ki  =  SEF;  p  =  pressure,  <X>  =  Elliptic  Integral  of  Second  Kind. 
See  Table  1  slide. 
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Model  3 
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Model  6b 

Center  Slice  (t  =  4.29  mm) 


Psf:  2.3x1  O'2  MPa 

Cf.-  175.30  mm 

af:  19.6  mm 

Data  zone:  (rave)2  -  (rave),  =  4.2635  -  0.4564  =  3 .807  mm 


Figure  4:  Unmultipli^;^  fi^mra:^||^^jfGe;#omiie  center 

of  the  cra(^fbivMb^er6  using. a  diffused  light  polariscope  for  a 


zone. 
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Mode  I  Algorithm  For  Determination  of  Stress  Intensity  Factor  (SIF) 

In  linear  elastic  fracture  mechanics  (LEFM)  using  the  photoelastic  approach,  one  can  begin  with  Mode  I  near  tip 
equations  (Fig.  A-l) 

~  +  =  n’2)  (A-l) 

where  Kj  is  the  Mode  I  SIF,  afj  are' the  contribution  of  the  nonsingular  stresses  in  the  measurement  zone,  and  r,6 
then  are  centered  at  the  crack  tip.  The  following  expression  is  computed,  in  truncated  form,  along  6  =  vr/2,  where 
fringe  spreading  is  greatest.  (Fig.  A-2)  Thus 


r 


nz 

max 


Kap  Kj 


+  r0 


(A-2) 


where  KAP  is  an  “apparent”  SIF,  which  includes  the  effect  of  a?-  {i.e.,  r0  =  /(a?-)}  with  the  singular  effect  in  the 

measurement  zone.  The  stress-optic  law  states  that  ,  where  N  is  the  measured  stress  fringe  order,  /  is  the 

material  fringe  value  and  t  the  slice  thickness.  Thus  T^zax  is  proportional  to  IV  and  may  be  regarded  as  the  measured 
quantity  together  with  r.  By  rearranging  terms  in  Eq.  A-2  and  normalizing,  we  can  obtain 


Kap$ 

Pv^ra  p-JWa 


(A-3) 


for  a  semi-elliptic  crack  where  the  coefficient  of  \J r j a  is  a  constant,  p  is  the  internal  pressure  and  a  is  the  crack 
size.  $  is  an  elliptic  integral  which  varies  with  the  aspect  ratio  of  the  crack  (a/c).  Its  form  is  approximated  by  yfQ 
where  Q  is  given  in  Table  I.  In  general,  when  applied  to  cylindrical  vessels,  the  denominator  of  Eq.  A-3  should  be 
P: t-  However,  in  the  present  problem  geometry,  Rf/t  =  1,  and  the  R{/t  can  be  dropped  here. 

By  defining  the  normalized  SIF  as 


p  _  KapVQ 
p<y¥a 


(A-4) 


one  can  plot  F  vs.  -JrJa  and  locate  the  linear  zone  implied  by  Eq.  A-3,  which  is  the  zone  dominated  by  the  stress 
singularity.  By  extending  this  line  to  the  origin,  the  value  of  F  is  determined  as  shown  in  Fig.  A.3  for  Model  1. 


Fig.  A-l:  Mode  I  Near  Tip  Notation  Flg'  A ~2''  M°de  1  Fhnge  (Pattern  Unmultiplied) 
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a/c  =  0.6  +  0.02 
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Figure  A-3:  Determination  of  Normalized  Stress  Intensity  Factor  (F)  from 
Test  Data. 


These  crack  shapes  were  not  semi-elliptical  or  through  the  length  cracks. 
Cylinder  length  was  336mm 


Summary 
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